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ABSTRACT 

We present first results of a study of the submillimetre (rest frame far-infrared) properties 
of z ~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galax- 
ies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories 
Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of 
LBGs well below the current limit of the survey. Although LBGs are not detected individ- 
ually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 /im) yields a 
statistically significant submm detection with mean flux (£250)= 5.9±1.4 mJy confirming the 
power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. 
In comparison, the Spitzer 24 fim detected BX/BM galaxies appear fainter with a stacked 
value of (S250) = 2.7±0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we de- 
rive median infrared luminosities, Ltr, of 2.8 x 10 12 L Q and 1.5 x 1O 11 L for z ^3 LBGs and 
BX/BMs, respectively. We find that Ltr estimates derived from present measurements are in 
good agreement with those based on UV data for z^2 BX/BM galaxies, unlike the case for 
z ~3 infrared luminous LBGs where the UV underestimates the true Lir. Although sample 
selection effects may influence this result we suggest that differences in physical properties 
(such as morphologies, dust distribution and extent of star-forming regions) between z ~3 
LBGs and z^2 BX/BMs may also play a significant role. 

Key words: Galaxies: high redshift -galaxies: starburst -submillimetre 



1 INTRODUCTION 

The broadband dropout technique has been a very succesful tool 
for discovering high-redshift galaxies, the so-called Lyman break 
galaxies (LBGs, e.g. Steidel & Hamilton 1993, Steidel et al. 1999). 
The initial selection focused on z ~3 samples. The same colour 
criteria were later extended to select Lyman break galaxies at 
1.4<z<2.5 (the so called BX/BM objects) with approximately the 
same range of UV luminosity and intrinsic UV colours as the z ~3 
LBGs (Reddy et al. 2006). The dropout broadband technique pro- 
vides a complete census of UV light at high redshift, with well 
over a thousand galaxies detected at z>1.5. Recent detailed stud- 
ies including Spitzer observations have shown that some of these 
galaxies have large stellar masses >10 10 Mq (e.g. Rigopoulou et 
al. 2006, Reddy et al. 2006, Magdis et al. 2008, 2010a) while their 
comoving volume density at z ~3 is ~0.005 Mpc -3 (e.g. Reddy 
& Steidel 2009). 

A number of issues related to the nature and properties of 
z ~3 LBGs remain unclear. The dust-corrected star formation rate 
(SFR) of LBGs can be as high as 100 M0/yr, which would cor- 
respond to Ssso ~ lmJy depending on specific dust parameters 
(Chapman et al. 2009). However, the search for the sub-millimetre 
(submm) counterparts of LBGs has proven challenging due to un- 
certainties in the relations used to predict the rest-frame far-infrared 
luminosity from the UV. Peacock et al. (2000) analysed the submm 
emission from star forming galaxies with the highest UV star- 
formation rates and found that they were statistically detected with 
a flux density Ss50=0.2 mJy for a star formation rate of 1 h -2 
MQ/yr. Chapman et al. (2000, 2009) reported the submm detec- 
tion of Westphal MMD-1 1 and Westphal-MM8, while Rigopoulou 
et al. (2010) reported mm detections of a further two LBGs, EGS- 
D49 and EGS-M28 selected based on their strong MIPS 24 /im 
emission (e.g. Huang et al. 2005). Despite these promising detec- 
tions the properties of the FIR and submm emission from LBGs, 
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their dust content and their possible contribution to the cosmic far- 
infrared background is still largely unconstrained. 

With the advent of Herschel (Pilbratt et al. 2010) it is now 
possible to investigate the submm (rest-frame far-infrared) proper- 
ties of LBGs. In this letter we report first results on the far-infrared 
properties of LBGs based on observations that are part of the Her- 
schel Multi-tiered Extragalactic Survey (HerMES), a Guaranteed 
Time project that will eventually result in a variety of surveys of 
varying depth and area which will be covered in five photometric 
bands (110, 160, 250, 350, 500 /Ltm; Oliver et al. 2010). The re- 
sults presented here are based on HerMES data taken as part of the 
Herschel Science Demonstration Phase. Throughout this paper we 
assume S7 m = 0.3, 57a=0.72 and Ho = 72 km s -1 Mpc -1 . 



2 OBSERVATIONS, SAMPLE SELECTION AND 
ANALYSIS 

2.1 Herschel Observations 

Submm observations of the Northern field of the Great Observa- 
tories Origins Deep Survey (GOODS-N) were carried out at 250, 
350 and 500 /im, with the Spectral and Photometric Imaging Re- 
ceiver (SPIRE). The instrument and its capabilities are described 
in Griffin et al. (2010), while the SPIRE astronomical calibration 
methods and accuracy are outlined in Swinyard et al. (2010). The 
GOODS-N images are amongst the deepest possible with SPIRE 
and, the instrumental noise is less than the confusion noise from 
overlapping faint sources. Confusion noise values of 5.8, 6.3 and 
6.8 mly beam -1 at 250, 350 and 500 ^im respectively, are reported 
in Ngyuen et al. (2010). Besides blind source extraction resulting 
in single-band catalogues (SCAT, see Smith et al. 2010, in prep.), a 
novel source extraction method based on 24 /im priors has been de- 
veloped to detect sources as close as possible to the confusion limit 
(see Roseboom et al. 2010, hereafter XID catalogue). The method 
uses a matrix inversion technique which relies on the assumption 
that sources detected in the 250 /im band will also be detected at 
24 fim deep surveys (e.g. Marsden et al. 2009). The 24 /im cata- 
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logue positions are then used to find sources in the Herschel 250 
(im images. The flux densities of the sources are allowed to vary 
until finally a set of flux densities is found that produce the best 
match to the image. In the current study we have made use of both 
blind (SCAT) and 24 /im prior source catalogues (XID), while for 
the stacking analysis we have used calibrated GOODS-N SPIRE 
images. 

2.2 Sample Selection and Analysis 

The GOODS-N region contains 58 UV selected z ~3 LBGs (Stei- 
del et al. 2003) and 212 UV selected BX/BM objects (Reddy et al. 
2006) down to R s^25.5. 49 LBGs and 200 BX/BMs have been 
detected with the Spitzer Infrared Array Camera (IRAC, (down 
to 25.0 mag(AB) at 3.6 fira) and 9 LBGs and 69 BX/BMs have 
also been detected with the Multi Imager Photometer for Spitzer 
(MIPS, down to S24=20 /iJy, 5a). In the current study we investi- 
gate the submm properties of the UV selected LBGs and BX/BMs 
focusing particularly on the sub-sample with MIPS detections: 9 
z -3 LBGs and 69 1.5< z ^2.5 BX/BM (hereafter the MIPS- 
LBG and MIPS-BX/BM samples). All LBG and BX/BM galaxies 
have spectroscopic redshifts determined from optical spectroscopy 
(Steidel et al. 2003, Reddy et al. 2006) which has also been used 
to confirm the absence of strong high ionization emission lines in- 
dicative of the presence of AGN. Objects classified as AGN/QSO 
are excluded from this study. 

The LBG and BX/BM samples were matched to the blind 
(SCAT) and priors (XID) catalogues. An object is considered de- 
tected when its flux is at least 3a above the noise (confusion plus 
instrumental). None of the LBGs and only three of the BX/BM ob- 
jects are detected in the HerMES GOODS-N priors catalogue down 
to S250 ~20 mjy. We discuss the properties of individual sources 
in section 2.3. 

To assess the reliability of SPIRE detections at faint flux levels 
«20 mjy) close to the confusion limit we first look at the flux dis- 
tribution (per pixel) of the whole GOODS-N map. For the present 
discussion we focus on the 250 /im band which is the most sensi- 
tive and has the smallest beamsize (18")- We find that the fraction 
of pixels above 5, 10 and 20 mJy is 6, 4 and 0.6 percent, respec- 
tively. These percentages imply that 16(4.5), 10(3) and 1.5(0.23) 
sources out of the 270(78) sources studied here could be associ- 
ated with spurious detections at 5, 10 and 20 mjy. The number in 
parentheses correspond to the MIPS detected objects. 

Since none of the LBGs and only three of the BX/BMs are 
detected individually, we examine the properties of the two samples 
via stacking analysis. We measure flux densities directly from the 
calibrated GOODS-N SPIRE 250, 350 and 500 /im images (at the 
optical positions of the LBGs and BX/BMs). For the stacking we 
first consider the optically-selected LBG and BX/BM samples. For 
the analysis we employ median stacking and stack at the optical 
position for each object. We exclude sources near bright objects 
(within 18" at 250 /im) to avoid contamination of the measured 
signal. Stacking the UV selected samples returned no detections in 
any of the three SPIRE bands. For the LBGs we determine 3a upper 
limits of S250 <2.8 mjy, S350 <1.5mJy and S500 <0.9mJy. For 
the BX/BMs the corresponding (3a) upper limits are S250 <2.02 
mjy, 5350 <1.2mJy and S500 <0.6mJy, respectively. 

The same stacking technique was employed to investigate the 
properties of the MIPS detected LBGs and BX/BMs. In Figure 
1 we show the histograms of the 250 /im flux density distribu- 
tions for the MIPS-LBGS and MIPS-BX/BMs. In the same plot, 
we show the normalised distribution of flux densities per pixel for 



the whole GOODS-N image. In both cases the distributions have 
positive skews and in the case of LBGs there is a clear positive tail 
implying that pixels containing flux from the LBGs have higher 
flux than the average pixel in the map. The case for the BX/BMs 
is not as clear. Although the BX/BM sample is bigger (69 objects 
vs 9) the population as a whole does not appear to strongly emit 
in the submm. In order to confirm that the two distributions are in 
fact different, we carry out a K-S two-sample test. The test results, 
D= 0.31 and (probability) a(D) = 0.03, suggest that the two flux 
density distributions are intrinsically different at the 2a level. The 
mean flux densities are (S250) = 5.9±1.4 mjy (LBG) and {S250) 
=2.7±0.8 mJy (BX/BM). The errors reported in these measure- 
ments have been quantified by stacking at 9 (for the LBGs) and 69 
(for the BX/BMs) random positions and then repeating the process 
40,000 times. The la value of the derived distribution is adopted as 
the uncertainty of the measurement. Stacking at 350 and 500 /im 
resulted in no formal detections. Instead we quote 3a upper lim- 
its of S350 < 4.9 mjy and S 500 < 3.4 mjy for the MIPS-LBGs 
and S350 < 2.6 mjy and S 50 o < 1.6mJy for the MIPS-BX/BMs, 
respectively. 

A likely physical explanation for this difference in detection 
rates comes from the MIR properties of the two samples. For z ~3 
LBGs the MIPS 24/im band contains contributions from 6/im (hot) 
dust continuum plus the 6.2 /im emission from Polycyclic Aro- 
matic Hydrocarbons (PAHs). The PAH contribution to the MIPS 
24/im flux is ~ 20-30 percent (e.g. Huang et al. 2005), therefore, 
the MIPS 24 /im flux mostly reflects the strength of the underlying 
continuum. Thus, the MIPS detected LBGs (with the strongest 6 
/im continuum) are likely to be amongst the most luminous z ~3 
LBGs. For z~2 BX/BMs the MIPS 24/im band includes contri- 
butions from both the 7.7/im PAH feature and the underlying con- 
tinuum. However, z~2 galaxy populations show a wide variety in 
their L^s.s^m rest-frame luminosity, with BX/BMs, in particular 
showing a relatively weak MIR continuum (Reddy 2006b). PAHs in 
BX/BMs, might be instrinsically weak since metallicity is known 
to affect their strength (Houck et al. 2005). The combination of low 
MIR continuum plus weaker PAH features may be responsible for 
the low detection rate in the submm of 24 /im selected BX/BMs. 
We thus conclude that, z ~3 LBGs detected by MIPS are likely to 
be on average more luminous than MIPS detected z ~2 BX/BMs. 



2.3 Individual detections 

In the previous section we examined the average submm properties 
of LBGs and BX/BMs here we take a closer look at individual 
detections, focusing on objects with MIPS detections. As discussed 
in section 2.1 of the 69 MIPS-BX/BM objects, 3 are detected with 
S250 > 20mJy and S/N > 3 however, two of those, BX1296 and 
BX1223 lie close (within 3") to bright submm galaxies, BX1296 
to GN07 and BX1223 to GN06 (the notation is from Pope et al. 
2005). These confused cases have not been considered in this work. 
BM1326 is clearly detected, with 5 , 2so : =22±5 mjy while a further 3 
BX/BMs appear in the 250 SPIRE map (and prior-based catalogue) 
but with 250 /im flux densities < 10 mjy. Additionally, 2 (of the 9) 
LBGs, HDFN-M18 and HDFN-M23 appear to be present in the 
maps (and prior-based catalogue) although with fluxes below the 
lOmJy level. We note that, HDFN-M23 is included in the 5a radio 
catalogue of Morrison et al. (2010), with a flux density 21.2±4 /iJy. 
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Figure 1. Histograms showing the 250 (im flux density distributions for LBGs (left) and BX/BMs (right). The blue curves show the flux density distribution 
in the entire GOODS-N field (control sample). In each case the control sample has been normalised by a factor equal to the total number of objects in each 
case (N=9 for LBGs, N=69 for BX/BM) divided by the total number of pixels in the whole GOODS-N image. 
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Figure 2. Rest frame average SED of MlPS-detected LBGs (left) and 1< z <2.5 BX/BMs (right) galaxies. For the SEDs we used mean values of UGRVUK, 
IRAC, MIPS, the mean value derived from SPIRE 250/jm and upper limits from SPIRE 350, 500 fim and Aztec 1.1 mm measurements (for z ~3 LBGs only). 
The rest-frame UV/optical is fit with BC03 models (magenta) while the mid/far-infrared part of the SED is fit with CE01 templates (green line). 



3 RESULTS 

3.1 Spectral Energy Distributions of LBGs 

Fig. 2 shows the average Spectral Energy Distribution (SED) of 
MlPS-detected LBGs and BX/BM galaxies. The SEDs have been 
constructed using available 'averaged' UGRVUK, IRAC, MIPS and 
SPIRE flux measurements. For the MIPS-LBGs we also use the 
1.1 mm Aztec measurement from Magdis et al (2010b). We fit 
the optical/near-infrared part with model SEDs generated using 
the Bruzual & Chariot (2003, BC03) code, while the mid-to-far 
infrared part is fit using Chary & Elbaz (2001, CE01) template 
SEDs. In brief, we use BC03 and construct stellar population mod- 
els with a Salpeter IMF and constant star formation rate, which has 
been shown (van Dokkum et al. 2004, Rigopoulou et al. 2006, Lai 
et al. 2007) to provide an adequate description of the properties 
of high redshift galaxies with ongoing star formation. Age, stel- 
lar mass, dust reddening E(B-V) and star formation rates are then 



derived from the model fits. It is beyond the scope of the present 
work to discuss these results, a detailed analysis of the properties 
of the stellar population in Spitzer detected LBGs can be found in 
e.g. Rigopoulou et al. (2006), Magdis et al. (2010a, for LBGs) and 
Reddy et al. (2006, for BX/BMs). It is however, worth noting the 
differences in the optical part of the SED with the BX/BM galaxies 
showing a much 'bluer' SED. 

We fit the far-IR/submm part with templates from the CE01 
library, with the best-fit templates rendering mean (I/ir) values of 
2.8(±0.6)xl0 12 L Q for MIPS-LBGs and 1.5(±0.5)x 10 11 L for 
MIPS-BX/BMs. The derived averaged L m for LBGs is typical of 
those seen in Ultraluminous Infrared Galaxies (ULIRGs). Using 
the infrared luminosities we derive average Star Formation Rates 
(SFRs) of 296M /yr and 245M /yr, for the z ~3 LBGs and 
BX/BM galaxies, respectively. The <SFR> derived from the IR 
for LBGs is in agreement with the radio SFR estimate (280±85Mq 
yr" 1 ) but higher than the UV SFR estimate (250+go M yr -1 ) 
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Figure 3. Modified blackbody temperature fit for the galaxy BM1326 
(which is individually detected at the SPIRE bands). The black filled 
squares represent the individual detections at 250, 350 and 500 fim. For 
the fit we have fixed /3 c j=l .5. 
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Figure 4. A comparison of estimates of Lir from the present submm obser- 
vations and from the UV (Lib, UV f° r the average z ~3 LBGs (magenta), 
the average z ~2 BX/BMs (green), together with the BM1326 and D49 
and M30 (the two LBGs with 1.2 mm detections). The solid line represents 
the Lib.,uv ~ -km- 



reported in Magdis et al (2010b). Turning to the BX/BM galaxies, 
the present SFR estimates agree well with those derived from the 
UV (Reddy et al. 2006) for UGR-selected galaxies. 



3.2 Dust, temperature and mass 

To derive the dust temperature, we use a single temperature grey- 
body fitting function (Hildebrand 1983) in which the thermal dust 
spectrum is approximated by: F v — Q V B V (T&), where B v is 
the Planck function, Q v — Qo(i / /i / o) l3d is the dust emissivity 
(with 1$: fid ^2) and, Td is the effective dust temperature. For 
hu I kTd ^ 1 the spectrum becomes : 



F v oc 



exp{hv / kTd) 



(1) 



Since Td and fid are degenerate for sparsely sampled SEDs we have 
fixed fid = 1.5 (e.g. Blain et al. 2003) which is consistent with 
SED fitting of low and high-z systems (e.g. Dunne et al. 2001). 
A higher value of fid will result in lower dust temperatures (Sajina 
et al. 2006). The dust temperature for BM1326, was obtained from 
the best fit model derived from minimization of the \ 2 values. The 
uncertainty in the measurement was obtained by repeating the pro- 
cedure based on perturbed values of the photometric points within 
their errors. To derive dust masses we follow: 



S v Dl 



K(\rest)Bv(\ res t, Td) 



(2) 



where Aid is the total dust mass, S„ the observed flux density, 
Dl luminosity distance, n(X r est) is the rest frame dust mass ab- 
sorption coefficient (taken from Weingartner & Draine 2001) and 
B 1 /(A reS i, Td) is the Planck function. For the z ~3 LBGs we 
assume Td = 45 K, a value chosen from Td estimates for lo- 
cal ULIRGs (Lisenfeld, Isaak, Hills et al. 1997) since the aver- 
age MIPS-LBG appears to have Lir > 10 12 Lq. We derive dust 
masses of A/ d =5.5±1.6x 10 s M and ALj=12.8±2.3x 10 s M Q for 
the LBGs and for BM1326, respectively. 



4 DUST OBSCURATION IN UV SELECTED GALAXIES 

The present SPIRE observations allow us to probe the cold dust 
peak of LBGs, determine their far-IR luminosity, dust temperature 
and dust mass from the far-IR alone with minimal additional as- 
sumptions. Earlier attempts to detect submm emission from LBGs 
(with targets selected mostly based on SFR(UV) estimates) were 
not met with success (e.g. Chapman et al. 2000, Peacock et al. 
2000). These initial results led to suggestions that either Td is high 
(Td ^90 K) or, that estimates of Lir from the rest-frame UV 
and/or from the scatter in the UV-slope/far-IR relation are uncer- 
tain (e.g. Chapman et al 2000). Recently, Rigopoulou et al. (2010) 
reported 1.2 mm detections of 2 LBGs in the Extended Groth Strip 
(EGS), both are detected in the MIPS 24 /im imaging survey of the 
EGS (for the full SEDs see Rigopoulou et al 2006). Briefly, their 
properties are similar to those of the GOODS-N z ~3 LBGs with 
5*24 in the 50-100 /iJy range. Using CE01 models we infer infrared 
luminosities, Lir ~ fewx 10 12 Lq for each of these LBGs. 

Let us now focus on the MlPS-detected LBGs and BX/BMs 
using the mean LBG properties reported in section 3.1 and com- 
pare Lir values from the SPIRE data to those derived from the UV. 
Lir,uv is determined as follows: at z ~3, G and R bands corre- 
spond to rest-frame 1200 A and 1500 A, respectively, thus alow- 
ing us to estimate the slope fi. Assuming solar metallicity, Salpeter 
IMF and continuous dust-free star formation models we use the 
BC03 code to generate SEDs to fit each of the LBGs, assuming the 
Calzetti (2000) attenuation law (but see also Buat et al. 2010 for 
a discussion of alternative extinction laws). Based on the best-fit 
model we derive extinction values E(B-V) and infer the observed 
and intrinsic 1500 A flux density and subsequently, L1500 lumi- 
nosity. We repeat the same process for the two galaxies with mm 
detections and the z-~2 mean BX/BM galaxies (using the B band 
flux density to estimate the intrinsic L1500 luminosity). 

The results are plotted in Fig. 4. The UV appears to underes- 
timate the Lm : uv of both the averaged z ~3 LBG (and the two 
LBGs with additional mm detections) by a factor ~2. This is per- 
haps not surprising given that the LBGs that appear to be detected 
in the submm regime are all ULIRGs. It is known that the UV un- 
derpredicts the Lir for both local ULIRGs (e.g. Howell et al. 2010) 
and z>2 submm-luminous galaxies. On the other hand, it appears 
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Figure 5. Bolometric luminosity, approximated as the sum of the IR and UV 
luminosities, vs. IR-to-UV luminosity ratio (dust obscuration). Small blue 
circles are z ~2 spectroscopically confirmed BX/BMs (from Reddy et al. 
2006), red triangles represent submm-luminous galaxies, black crosses are 
local normal (Bell et al. 2003) and starbursts (Brandl et al. 2006) galaxies 
and cyan stars are ULIRGs. The large black circle corresponds to BM1326 
while the large green circle is the average BX/BM. Magenta and black 
squares are the average LBGs and the two LBGs with mm detections. The 
solid line indicates the best-fit linear relation for spectroscopically con- 
firmed (/GS-galaxies detected at 24 fj,m (from Reddy et al. 2006). The red 
and blue dashed lines are lines of constant Lijv and Lu luminosity. The 
errorbars for the stacked LBG (magenta) and BX/BM (green) values have 
been magnified for clarity. 

that the UV provides a better estimate (closer to the measured Lm) 
for the averaged BX/BMs. 

Finally, it is instructive to look at variations of the obscura- 
tion of these UV-selected galaxies. For this purpose we examine 
the bolometric luminosity (defined as the sum of the IR and UV 
luminosities) as a function of obscuration (approximated by the ra- 
tio of IR-to-UV luminosity) for LBGs and BX/BMs. For compari- 
son we include submm-luminous and UGR-selected z ~2 galaxies 
(from Reddy et al. 2006 and references therein). The resulting plot 
is shown in Fig. 5. The straight line indicates the correlation found 
by Reddy et al. (2006, 2010) for z ~2 UGR-selected galaxies. The 
averaged z ~3 LBG and the two individually detected ones appear 
to follow the relation defined for the z~2 galaxies. In terms of lumi- 
nosities, both averaged LBGs and BX/BMs have similar Luv (few 
x \Q W Lq) but LBGs have higher Lm and thus higher Lfir/ Luv 
ratio. Since it is well established that obscuration decreases with in- 
creasing redshift (Reddy et al. 2006, 2010, Adelberger and Steidel 
2000), the difference in the Lfir/ Luv ratio must be attributed to 
different causes. While selection effects are likely to play a role (see 
section 2.2) we argue that possible differences in morphologies, 
dust distribution and extent of star-forming regions are also likely 
to contribute. Morphological studies of UV-selected z ~2 and z~3 
galaxies in the GOODS-N field find few differences between the 
two samples (Law et al. 2007) although, dustier galaxies (as evi- 
denced by E(B-V)) were found to show more nebulous UV mor- 
phologies. Finally, since MlPS-detected LBGs have ULIRG-like 
luminosities (section 3.1) it is possible that their UV and IR emis- 
sion originates in different regions (as observed in local ULIRGs 
e.g. Wang et al 2004, see also Huang et al. 2007) and thus could 
account for the higher Lir/Luv ratio observed. 
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